www.gusucode.com > 调制识别 基于决策论的通信信号调制识别实现matlab源码程序 > moulation classification/key feature simulink/M1_fig.m
clear;fb=500;fc=2000;fs=12000;Ac=1;f1=500;f2=1000;f3=1500;f4=2000; for i=1:30 for k=1:100 k d2=randint(1,300,2);d4=randint(1,300,4);d16=randint(1,300,16);d8=randint(1,300,8); s1=ask2(d2,fb,fc,fs,Ac);s2=ask4(d4,fb,fc,fs,Ac);s3=qam16(d16,fb,fc,fs,Ac); s4=fsk4(d4,fb,f1,f2,f3,f4,fs,Ac);s5=psk8(d8,fb,fc,fs,Ac); ss1=awgn(s1,i,'measured');ss2=awgn(s2,i,'measured'); ss3=awgn(s3,i,'measured');ss4=awgn(s4,i,'measured'); ss5=awgn(s5,i,'measured'); Ma11(k)=M1_func(ss1);Ma12(k)=M1_func(ss2);Ma13(k)=M1_func(ss3); Ma14(k)=M1_func(ss4);Ma15(k)=M1_func(ss5); end M11(i)=mean(Ma11);M12(i)=mean(Ma12);M13(i)=mean(Ma13); M14(i)=mean(Ma14);M15(i)=mean(Ma15); end plot(M11,'-*') hold on plot(M12,'-^') plot(M13,'-s') plot(M14,'-d') plot(M15,'-p') xlabel('snr');ylabel('M1'); grid legend('ask2','ask4','qam16','fsk4','psk4',0);